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Abstract The cosmic microwave background (CMB) represents a unique source for 
the study of gravitational lensing. It is extended across the entire sky, partially polar- 
ized, located at the extreme distance of z = 1 100, and is thought to have the simple, 
underlying statistics of a Gaussian random field. Here we review the weak lensing 
of the CMB, highlighting the aspects which differentiate it from the weak lensing 
of other sources, such as galaxies. We discuss the statistics of the lensing deflection 
field which remaps the CMB, and the corresponding effect on the power spectra. We 
then focus on methods for reconstructing the lensing deflections, describing efficient 
quadratic maximum-likelihood estimators and delensing. We end by reviewing recent 
detections and observational prospects. 
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1 Introduction 

1.1 What is the CMB ? 

The cosmic microwave background (CMB) is a blackbody radiation which permeates 
the Universe. It was released approximately 400,000 years after the Big Bang, when 
the primordial plasma had cooled enough that neutral hydrogen atoms could form, 
and the Universe became effectively transparent to radiation. To a very good degree 
of approximation, the CMB temperature is uniform across the sky, with a value today 
of 2.725 K [[!]]. At the level of one part in 10 5 , however, the CMB temperature has 
small anisotropics which are due to fluctuations of density and bulk velocity in the 
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primordial plasma. The measurement of these anisotropies is clearly a difficult task, 
but a number of very successful experiments have driven rapid progress in this field 
over the past two decades. Modern temperature experiments such as ACtQ, SPT^j, and 
Planck^ are now pushing to robust arcminute-scale measurements, while polarization 
experiments like Spider^, QUIET^j, BICEP2 and Planck[] will soon obtain low-noise 
maps of the CMB polarization. Gravitational lensing effects should become readily 
apparent in both types of observation, requiring detailed modelling to interpret the 
data correctly, but also offering additional information about the Universe at interme- 
diate redshifts. 



1.2 What is special about the CMB in the context of lensing? 

In the context of lensing, the CMB can be thought of as a single source "plane" with 
several unique features: 

1 . It is extended across the entire sky. 

2. It probes the matter distribution at relatively higher redshifts (z ~ 2) than can be 
studied with e.g. galaxy lensing: the last-scattering surface from which the CMB 
was released (z ~ 1100) is the most distant light source we will ever be able to 
see. 

3. It is partially linearly polarized, which provides two additional lensed observ- 
ables. The local quadrupole of the radiation field at the last scattering surface 
generates linear polarization through the mechanism of Thomson scattering, at a 
level of 10% of the temperature anisotropies [g]. The observable CMB is there- 
fore characterized by three numbers as a function of position on the sky - the 
total radiation intensity, and two Stokes parameters describing the amplitude and 
direction of the linear polarization. 

4. To a good degree of approximation the source fluctuations may be treated as those 
of a Gaussian random field. The signal along any particular direction is essentially 
random, but the correlations between directions at different angular separations 
is encoded in the power spectrum, which completely characterizes the signal. A 
typical CMB temperature realization, for example, is plotted in Fig. |l| 



1.3 Why study CMB lensing? 

To someone interested in the primary CMB, the lensing of the CMB is important as 
a contaminant. As we will discuss in Sec. lensing has a small but non-negligible 
effect on the CMB power spectra, and must be accounted for in modern parameter 

1 ^tt^jJ_/j^]^_ : J2^l^2^-L212^£^ on • edu/act/ 

2 itt^j^/gole^ichicago^edu/ 

3 ^^p^^W^W^^^C^^^^^^^^^texjDhj3?jDr^^e^t=p^^a]^l* 

4 ittp : //www, astro . caltech.e du/- lgg/spider_f ront .htm 

5 ittp : //quiet . uchicago . edu/ 

6 ittp : //www. rssd. esa. int/ index. php?project=plancl? 
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Fig. 1 Left: lensed CMB realization. Middle: difference map between lensed and unlensed CMB. Right: 
realization of \y for the lensing, and an overlay of its gradient, the deflection angle. 



analyses to obtain unbiased constraints. Perhaps more importantly, lensing effects 
generate a curl-like (B mode) polarization pattern on the sky which acts as a limiting 
source of confusion for low-noise polarization experiments targeting the signal from 
primordial gravitational waves This confusion can be reduced with an accurate 
cleaning of the lensing-induced signal, which we will discuss in Sec. |5J 

Apart from being a nuisance for traditional observables, the lensing of the CMB 
can act as an additional source of information. A typical analysis of the CMB assumes 
Gaussianity and statistical isotropy, in which case the power spectrum is the only 
quantity of interest. As we shall discuss, lensing can be thought of as introducing into 
the CMB small amounts of non-Gaussianity (when marginalized over realizations of 
the lenses) or statistical anisotropy (for a fixed distribution of lenses). This effectively 
introduces information into the CMB, contained in the higher-order statistics (for 
the non-Gaussian viewpoint) and the off-diagonal elements of its covariance matrix 
(for the anisotropy viewpoint). With only one CMB sky to observe and interpret, 
both these viewpoints are useful. The additional information from lensing probes the 
state of the Universe at intermediate redshifts (z ~ 2). This can be used to break 
parameter degeneracies and place improved constraints on quantities that affect the 
geometry or density perturbations at late times, such as the dark energy equation of 
state and portion of the energy budget in massive neutrinos. An optimal analysis of 
lensing effects with the data from the Planck satellite, for example, will enable us 
to measure the sum of neutrino masses to ~ 0.1 eV, while lens reconstruction with a 
next-generation polarization mission such as EPIC/CMBPol can constrain the sum to 
0.05 eV or better [||,[|[7]]. This is an interesting limit, close to the minimum value for 
the sum of the masses suggested by terrestrial oscillation measurements in the normal 
hierarchy. 
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2 Observables 

Here we will briefly review the mathematical description of CMB observables and 
establish our notation. Throughout this review we will use an overtilde to repre- 
sent lensed quantities For quantitative calculations we assume that the Universe is 
described by a flat, concordance ACDM cosmology with parameters given by the 
WMAP 5-year best-fit model. 

The CMB radiation pattern on the sky is conveniently expressed in terms of the 
Stokes parameters T, Q and U, which are measures of the total intensity (T; here ex- 
pressed as an equivalent blackbody temperature relative to the background, uniform 
CMB) and linear polarization (Q and U). Rotation of the measurement coordinates 
does not affect T, however the Q and U Stokes parameters mix according to 

[Q±iU}^=e ±2i ^[Q±iU] , (1) 

where ^ is the angle of the rotation[| Thus, Q ± ill is a spin ±2 quantity. A rotation of 
45 degrees converts Q — > U, and a 90-degree rotation results in a sign flip. For more 
information on this notation, see for example |Q| . 

The statistics of the CMB are often simplest to study as a function of scale rather 
than angular position. On a patch of the sky which is small enough that it may be 
treated as flat, the Fourier modes provide an appropriate basis. In this review, we will 
take the convention that 

TO) = J ^T(x)e-** 
[E(l) ± W(l)] = - J^[Q± iU]e^e-^. (2) 

Here ^\ is the angle between 1 and the jc-axis. The placement of the e^ 2 '^ 1 factor in 
the polarization transform rotates the Fourier transforms of the Stokes parameters 
to a basis adapted to the Fourier wavevector, resulting in E and B fields which are 
spin-0. This makes E and B independent of any measurement conventions, and so 
they are the natural basis to perform calculations in. The E/B-mode decomposition 
is also useful for studying the origin of fluctuations in the early Universe: B-modes 
are not generated by scalar perturbations, and so the detection of non-zero primordial 
Z?-mode power could in principle be a "smoking gun" for primordial tensor modes 
0. 

The primary, unlensed CMB fluctuations are thought to be well approximated as 
a statistically-isotropic Gaussian random field with zero mean. This means that in 
the Fourier domain, each mode is drawn from a (complex) Gaussian, with variance 
which depends only on the length of the wavevector, and not on its direction. The 
length-dependent part is known as the power spectrum C/. In the Fourier basis, for 
example we have 

(X(l)Z*(l'))=S(l-l')Cf z , ( 3 ) 



7 This is for Stokes parameters defined on a right-handed x-y basis with z along the propagation direc- 
tion, and for a left-handed rotation about z. 
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where X, Z can be T, E or B, and / = |1|. If parity is respected in the mean, Cf and 
Cf B must vanish, but note that there is a non-zero cross-correlation between the tem- 
perature and the E'-mode polarization. Because the modes are drawn from a Gaussian 
distribution, they are completely characterized by their two-point statistics. Higher 
order odd-point functions are zero, and the even-point functions can be computed 
from the power spectra using Wick's theorem. 

We now proceed to describe the large-scale structures responsible for the lensing 
of the CMB. 

3 Lensing potential 

3.1 Lensing deflection 

The effect of lensing on the CMB can be described in the flat-sky limit as a remapping 
of the primary (unlensed) CMB given by 

f(x) = T(x + a(x)) (4) 
[Q ± id] (x) = [Q ± iU) (x + a(x)) (5) 

where a(x) is a field of deflection vectors. 

The observed deflection vector is given by a sum of all the deflections by lenses 
between ourselves and the last scattering surface; mathematically this is given by an 
integral over the standard weak lensing equation along the photon path. In a flat uni- 
verse, and under the Born approximation (taking the integral along the unperturbed 
path) we have 

a(x) = -2 ( X * dx^t V±f (*x; 770 - %), (6) 
Jo X* 

where f is the Weyl potential, x* is the conformal distance to last scattering, and rjo 
is the conformal time today. Here, the gradient is taken transverse to the propagation 
direction. In Einstein gravity with vanishing anisotropic stresses, the Weyl potential 
is simply the Newtonian potential. More generally, in any metric theory of gravity, f 
is given in terms of the Newtonian-gauge potentials 0n an d Ws by W = (0n + V / n)/2: 
it is the scalar field which determines the Weyl tensor (non-local part of the Riemann 
tensor) in the presence of scalar perturbations. For further details and the simple 
generalization to a non-flat Universe see JTl|]. 

3.2 Order of magnitude 

It is useful to get a feel for the order of magnitude of the lensing effect on the CMB, 
and why it is important. We can think of the large scale structure of the Universe 
as a collection of potential wells, crudely approximated as those due to point-source 
masses. Each gives an observed deflection angle of 

|a ps (x)|=4f*^W(*x), (7) 



6 




where % is the conformal distance to the point source and f ps (jx) is the potential at 
closest approach. Thus nearby lenses have the largest contribution from the geometric 
factor (x* — x) /X*> though the linear potentials are actually somewhat smaller at low 
redshift due to the suppression of perturbation growth by the dark energy. During 
matter domination the geometric factor is given by (1 — %/%*) ~ [1/(1 +z)] 1//2 , which 
depends somewhat weakly on redshift, so we will ignore it for the rest of this order- 
of-magnitude calculation. 

The potentials are nearly scale-invariant on large scales, but become small on 
scales below the turnover in the matter power spectrum. We therefore need to consider 
lensing by perturbations from the largest observable scale ~ x* ~ 14000 Mpc down to 
~ 300 Mpc. The number of lenses along the line of sight also depends on scale: there 
is only one of the largest lenses, but ~ 14000/300 ~ 50 of the smallest ones, so the 
smaller lenses contribute most to the lensing deflection. The size of each deflection 
is set by the Newtonian potentials, which have variance ~ (2.7 x 10~ 5 ) 2 per unit 
log range in scale, so typical deflections have \a\ ~ 10~ 4 . Taking each lens to be 
independent, the total r.m.s. deflection should therefore be ~ 50 1 / 2 x 10~ 4 ~ 7 x 
10 -4 , or about 2 arcmin. 

The coherence scale of the deflection field is given by the angular extent of a 
typical 300 Mpc structure. At a distance x = 7000 Mpc (halfway to the last-scattering 
surface) this is (300/7000) - 2°. 

These order of magnitude estimates are borne out by more accurate calculations. 
The RMS deflection angle for the WMAP 5-year cosmology is 2.7 arcmin. The scale 
dependence (i.e. angular power spectrum) and cumulative redshift dependence of 
the deflection power are shown in Fig. ||. Nearby objects tend to contribute to the 
lensing potential on larger angular scales. The deflection power peaks in the region 
20 < I < 100, corresponding to a coherence length of a few degrees. 
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3.3 Lensing potential 

Rather than using the deflection field directly, weak lensing studies are often cast in 
terms of the magnification matrix A, 7 - defined by 

A,^8 i]+ 4- aj =( i - K -^ r r2 + w V (8) 

1 1 dOi 1 \-72-(0 \-K+Y\J 

where the quantities K, j\ + iji, and co are known as the convergence, (complex) 
shear, and field rotation respectively. In galaxy lensing the magnification matrix is 
useful because in many cases this can be approximated as being constant over the 
angular extent of the source. Shear then produces an additional elliptical distortion of 
the projected shapes of galaxies, while convergence changes their apparent size and 
density on the sky. For galaxy lensing the shear is the more useful observable since 
observing the convergence would require knowledge of the unlensed distribution of 
galaxies. In CMB lensing, the degree-scale acoustic features are large compared to 
the coherence size of the magnification matrix, so the latter cannot be assumed con- 
stant over typical CMB features. Moreover, the statistics of the source plane are well 
known (the CMB is nearly a Gaussian random field), and so the convergence is just 
as useful as shear. For these reasons CMB lensing studies usually work directly with 
the lensing deflection field as a function of position on the sky, rather than using the 
magnification matrix. 

From Eq. (|[) we know that the deflection field is a sum of projected gradient 
terms. Approximating the light path by the unperturbed line of sight, the deflection 
field can be expressed as an angular gradient. In detail, relating the projected and 
angular gradients at distance % via V ± f = (Vf) / %, we can define a lensing potential 
y/(x) so that 

oc(x)=Vy(x), (9) 

where the derivative is now with respect to position on the sky. The lensing potential 
can be calculated as the line-of-sight integral 



X* 



X* X 



V(x) = -2 dx[^-^)^{xx;n -x)- (10) 
Jo \ X*X / 

To what extent is the deflection angle approximated by Eq. (|l0|) valid? The major 
assumption was the Born approximation, which takes the integral of Eq. (10) along 
the unperturbed path. This will be a good approximation if the deflection angles are 
sufficiently small, or if all of the lensing occurs in a single thin plane. The dependence 
of the lensing potential on this assumption can be investigated quantitatively under 
the second-order Born approximation - taking the gravitational potential along the 
perturbed path given by the first-order result. The effect on the gradient component 
of a(x) is computed in [|l2l|l3|], and is expected to be negligible for any future exper- 
iment. A potentially more important aspect of the second-order Born approximation 
is that it introduces a non-zero curl component into the deflection field, however this 
is not expected to be significant for experimental noise levels achievable in planned 
future surveys (i.e. greater than 0.25 /iK-arcmin) The validity of the Born ap- 
proximation can also potentially be investigated by ray-tracing through numerical 
simulations [O], although this has not yet been done. 
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3.4 Lensing potential power spectrum 

The lensing potential depends on the geometry of the universe and the gravitational 
potential along the line of sight, and hence, like galaxy lensing, is a probe of cosmol- 
ogy. Assuming linear evolution from the primordial density fluctuations, we can take 
f in Fourier space to be given by 

f(k,i7) = 7V(£;77)^(k), (11) 

where Tip (£;?]) is the transfer function, and M{k) are the primordial curvature fluc- 
tuations, with 3D power spectrum £5^(k). Under this assumption, expanding y/(x) 
in harmonic space and taking the angular power spectrum gives [111]] 



x* 

dxTy{k;rio-X)ji(kX) 



o V X*X 



X* X 



(12) 



where ji(k%) is a spherical Bessel function. The deflection power spectrum, 1(1 + 
1 )Cj ' , for a typical model is plotted in Fig. [| It is customary to plot [/(/ + 1 )] 2 Cj v / (2k) 
because 

[l(l + l)] 2 Cj v , x2 
dlnl^ — J 1 1 ^(a) 2 , (13) 
2n 

the variance of the deflections. 

The assumption of linear theory and Gaussian primordial fluctuations implies 
that f is Gaussian as well. However, on small scales it is clear that the potentials are 
subject to a large degree of non-linear evolution. Fortunately, the size threshold for 
non-linearity is small enough that there are many such non-linear clumps on the way 
to the last-scattering surface. This Gaussianizes the projection of f to some extent, 
but the non-linearity does still appear as an enhanced power spectrum, C, , on small 
scales; see Fig. |[ Non-linearity may therefore be approximately accounted for by 
scaling the linear transfer functions according to 



/ canon-linear/ k. „\ 

r '"^H «S <i4) 

The non-linear power spectrum can be calculated from analytical fitting formulae 
such as those of HALOFIT [[To|; this approach is implemented in the CAMB code 
Jl7|]. The ultimate test of the Gaussian approximation is to propagate CMB photons 
through large numerical A^-body simulations. Several groups have developed methods 
for creating the concentric-shell mass maps necessary to perform the radial integral 
of Eq. (||) numerically Jl5[|l8| ,|T9|], and all find that the HALOFIT approach works well 
at the percent level for I < 2000. 



4 Lensing effects on CMB observables 

In Fig. [j] we show the difference between the lensed and unlensed sky for a real- 
ization of the CMB temperature. In this picture, the effect of lensing is to create a 
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non-Gaussian "boiling soup" pattern in the observed CMB. Large-scale lenses create 
coherent displacements across large patches of the sky, which appear as a modulation 
of the gradient of the unlensed CMB in the difference map. 

To study the properties of the lensed CMB analytically, it is often useful to use 
a Taylor expansion of the lensing displacements. With the CMB temperature, for 
example, we have 

f(x) = T(x + a(x)) 

= T(x) + a a V a T(x) + ^a a a b V a V h T(x) + --- (15) 

The first line is not strictly correct outside the fiat-sky approximation due to the cur- 
vature of the sky, however the Taylor expansion of the second line is, provided that 
the gradients are treated as covariant derivatives. This is also true for polarization 
provided that the polarization tensor is used in place of T [p0||. 

The expansion in Eq. ( p"5| ) is frequently useful to gain intuition into lensing ef- 
fects. Consider, for example, the first-order expansion written in harmonic space: 

[a fl V fl r](l) = [V>V a T](l) = -/ ^L'-(1-L)va(L)T(I-L). (16) 

Here we can see that a mode of y(L) with wavevector L couples the unlensed CMB 
at wavevector 1 L into the observed CMB at wavevector 1. For fixed lenses, this 
mode-coupling introduces off-diagonal components into the covariance matrix of the 
observed CMB, with a characteristic spacing of 81 = 50 given by the peak of the 
deflection angle power spectrum. Specifically, 

(f(li)f(l 2 ))cMB = ^y(li+l 2 )- [hCf* +l 2 Cf 2 T ] (h^h), (17) 

where the average is over an ensemble of CMB fluctuations for fixed lenses. This 
is an interesting feature; large- wavelength modes of the lensing potential introduce 
correlations between CMB modes on much smaller scales. In Sec. || we will see how 
a reconstruction of the lensing potential can be made by correlating a large number 
of small-scale modes. 

We note that the convergence of the lensing Taylor expansion is somewhat suspect 
- the typical RMS deflection angle of ~ 3 arcmin is comparable to the 10 arcmin 
scale length of modes at / = 1000, and so the expansion converges slowly at the 
map level. Nevertheless, when calculating ensemble-averaged quantities the lensing 
Taylor expansion is often quite accurate. In this case, it is relative displacements 
rather than the actual lensing effects which are important, and the Taylor expansion 



captures these reasonably well [21 1. However the expansion is not accurate enough 



for a percent-level calculation of the lensed power spectrum. 



4. 1 Power spectra 



The CMB power spectra completely characterize the primary fluctuations for Gaus- 
sian, statistically-isotropic models. Even though the lensed sky is slightly non-Gaussian, 
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the lensed power spectra are still useful for comparison with observations. Here we 
discuss the calculation of the lensed power spectra using the Taylor-series expansion 
described above. A more accurate non-perturbative result may be derived using corre- 
lation functions [^2[^,|^|, however the Taylor approach gives more intuitive results 
and is less mathematically involved. 

We will begin with temperature. For the purposes of compact notation, it is useful 
to write the lensing Taylor expansion schematically as 

f(l) = T{l) + 8T{l) + 8 2 T{l) + --- , (18) 

where the power of delta indicates the order in the lensing potential. The power spec- 
trum to lowest order in Cj v then follows from 

Q ff 5(1-1') « Cf T 8(l-l') + (8T(l)8T*(i)) + (T(l)8 2 T* (l')) + (8 2 T(l)T*(V)) . 

(19) 

In the flat-sky approximation, the gradient operations are easily taken in Fourier 
space: 8T(l) is given by Eq. while [ |25| l 



8 2 t(\) = - 1 -J ^ J ^i 1 -(i 1 +i 2 -i)i 1 -i 2 r(ii)^(i 2 )^(ii+i2-i), (20) 

from which the power spectrum may be derived. Assuming that both T(l) and 
have the statistics of an isotropic, Gaussian random field we find: 

(T{l)8 2 T* (I')) + (S 2 T(l)T* (I')) = -l 2 RVcf T 8(1-1') (21) 
(ST(l)Sr(l')) = 5(1-1') J ^[ L .(l-L)] 2 C,^ L |Cr(22) 
where R v is half of the total mean-squared deflection, 

1 i2\ _ J_ / ., qvin -7 



= ^(|Vv/| 2 ) = ^ J jl 4 Cj"" ~ 3 x 10- 7 . (23) 

Assembling these results together, we find the lensed temperature power spectrum to 
first-order in Cj 

Cf T = (l-l 2 RV)Cf T + J^[l'.(l-l')] 2 C p r l , l Cj; T . (24) 

The second term, which represents a convolution of the unlensed temperature power 
spectrum with the lensing potential power spectrum, almost cancels the R v term over 
the acoustic part of the spectrum, leaving only a small blurring of the acoustic peaks 
with a characteristic kernel width of 81 ~ 50. This is illustrated in Fig. |[ where we 
plot the lensed and unlensed spectra. 

The unlensed CMB has very little power on small scales, I > 3000, due to dif- 
fusion damping. Deep enough into the damping tail, any non-negligible contribution 
to Eq. ( p4{ ) from the unlensed temperature power spectrum will be from modes with 
/' -C /, and so the lensed spectrum may simplified as [25] 26|] 

r H 2 V t HI' (]'\ A C TT 

r TT ~ rvv I n' \\ 2 r TT ~ flrVV i I l ' ~ i 2 ^V¥pT n ^ 

c > ~ c < 7w (1 - 1)c '' ~ /c ' J-^^~ lc > R > (25) 
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where R T is defined analogously to RV in Eq. (|23|), and is half the mean-squared 
gradient of the unlensed CMB. This causes an increase in the lensed CMB power 
of ~ 10% at I > 2000, as we can see in Fig. ||. Qualitatively, there are therefore 
two lensing effects on the temperature power spectrum: large-scale lenses smooth the 
peaks of the observed power spectrum, while small-scale lenses introduce small-scale 
power in the damping tail by lensing of the large-scale temperature gradient. 

We now turn to the lensed polarization power spectra. For simplicity, we will 
assume that B(l) = 0, i.e. that there are no primordial Z?-modes. In the flat-sky limit, 
lensing remaps the Stokes parameters defined on a fixed basis in the same way as for 
the temperature. The calculation is therefore similar to the temperature case, with the 
lensed spectra to lowest order in Cj v given by [ p5| ] 

<** = (l-l 2 RV)Ci E + + J ^ [l'.(l-l')] 2 C | n, | C? £ cos[2(^-^ 1 )] (26) 
Cp = + / [l'.fl-0] 2 CjS|<?* cofm-b)] (27) 

Cf S = Jj0j2 [I'-Cl-l')] 2 ^,^ sin 2 [2(^-5)], (28) 

where R v was defined in Eq. (^3j). Lensing of purely £-modes therefore generates 
B-mode polarization, even when there is none initially. 

The lensing of the £-mode power spectrum, and the TE spectrum, is qualitatively 
similar to that of temperature. There is a convolution which blurs the features of the 
spectrum and shifts power into the damping tail. Quantitatively, the magnitude of the 
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Fig. 4 The unlensed ii-mode power spectrum ?(/+ \)Cf E /2ll (top; dashed), lensing deflection power 
spectrum [1(1 + XYfCj^ (dot-dashed), and the contribution to the large-scale lensed B-mode power spec- 
trum from each log / given by half their product (thin solid). The dotted line is the large-scale white spec- 
trum /(/ + \)Cf^ /2n given by Eq. (H), which can be compared to the full numerical result (thick solid). 
The lensed B-mode power spectrum has a close-to-white spectrum at / -C 1000. 



lensing effects is larger for the E-modes than for temperature, because the peaks in 
the unlensed power spectrum are rather sharper, leading to lar ger fractional changes 
in the lensed power of near the acoustic peaks []2^,^4jj25[]. 

Lensing of the CMB polarization causes power to leak from E- to B-modes, and 
this may be an important contaminant in the search for primordial B-modes on large 
angular scales. Since the unlensed E-mode spectrum peaks around / ~ 1000 (see 
Fig. |), for / < 1000 we can take |1'| > |1| in Eq. (|§) giving 

= JLy d -Lfcyfcf, E . (29) 

This is independent of I, corresponding to a white-noise spectrum of B-modes. In 
Fig. Q we can see that the lensed B-modes are nearly white for I < 100, with the 



level predicted by Eq. (29). For the fiducial model we have Cf B - 2 x 10~VK 2 , 



corresponding to an instrumental white-noise level of ~ 5 /iK-arcmin. The next gen- 
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eration of polarization surveys should achieve noise levels close to this and lensing 
will therefore become an important source of confusion for them. 

To assess the amount of additional cosmological information in the lensed spectra 
over the unlensed, consider attempting to constrain departures of C, from some 
fiducial spectrum: 

cr=(i+f,)cr\M.. oo) 

An eigenmode analysis of the Fisher matrix for this problem can be used to identify 
modes of /j which are constrained well by data. For a cosmic-variance limited exper- 
iment to / max = 2000, there are essentially two such modes [^7|. From the TT and 
EE power spectra, there is a mode which peaks at / = 100 and tracks the dominant 
lensing power. These are the modes of Cj v which smooth the observed power spec- 
tra. Including B-mode polarization, there is a second mode which is much broader, 
reflecting the fact that the large scale B-mode power is sourced by lenses on all scales. 
This mode peaks at / = 200, but has a broad tail extending past / = 1000. 

Increasing the Z max of the temperature measurement beyond / = 2000, the small- 
scale modes of the lensing potential would ideally also be constrained by the increase 
in the small-scale lensed power which they cause. However, in this regime extra- 
Galactic foregrounds and secondary effects will provide a source of confusion which 
is, as yet, not well understood. 

The amplitude of the well-constrained Cj modes is somewhat dependent on cos- 
mology, and so they may be used to place improved constraints on parameters. Some 
care is required with the B-mode power spectrum due to the highly non-Gaussian na- 
ture of the lens-induced B-modes leading to correlations between estimates of Cf B at 



different scales [^7,28 29]]. On the whole, the power spectra are a rather sub-optimal 



statistic for extracting information from CMB lensing, however. Improved constraints 
can be made with estimators which are designed specifically to detect the lensing sig- 
nal; these are discussed in Sec. ||. 



4.2 Bispectra 

The bispectrum, or three-point correlation function is given by 

<f(l 1 )f(l 2 )f(l3))=B(l 1 ,l 2 ,l 3 )5(l 1 +l 2 + l3). (31) 

It is the simplest, lowest-order statistic which can be used to test for a skewed non- 
Gaussianity of the CMB fluctuations and so it has received much attention (see [ p0[ 
for a review). Considering the effects to second order in the lensing potential, the 
lensed bispectrum is given in terms of the unlensed bispectrum by 



B(l 1 ,l2,l3)5(li+l 2 + l3)=B(li,l2,l3)5(li+l2 + l3)+((7 , ai)7 , (l 2 )[V V r-Vr](l 3 )) 

+ (r(i 1 )[Vv/-vr](i 2 )[v V /-vr](i 3 )) 

+ i(r(l 1 )r(l 2 )[V>vVv a rV 6 r](l 3 ))+2perms) .(32) 
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One effect of the lensing terms is to generate a smoothing of features in the primordial 
bispectrum analogous to the effect in the power spectra This effect is not 

significant enough to affect the analysis of CMB data in the near future however [Ej]]. 

A more important effect is due to the correlation of the CMB lensing potential 
with the large-scale temperature anisotropics. The CMB radiation we observe is red- 
shifted or blueshifted by the evolution of potentials along the line of sight, a phe- 
nomenon often referred to as the (late-time) integrated Sachs-Wolfe (ISW) effect $5% 
|33| ]. The ISW contribution to the CMB anisotropics from linear potentials is signifi- 
cant at low multipoles (£ < 100) due to the global effects of dark energy, but dies off 
rapidly on smaller scales. Since the evolving potentials that generate the ISW effect 
also cause lensing, there is some cross-correlation C Tv between T and which is 



calculated in e.g. [25|. The correlation can be neglected for calculating the lensed 



power spectra, however a non-zero cj^ generates a bispect rum, b y the second term 



in Eq. (|32[), even if the primordial fluctuations are Gaussian [26 34Q . At leading order 
in y/ this is given by 

(f(li)f(l 2 )f(l 3 )) «-^5(li+l 2 + l 3 )(C / : ' v 'C / 7l 1 -l2 + 5perms.). (33) 

The value of this induced non-Gaussianity is greatest for "squeezed" triangles in 
which (] ~ I2 > h ( or permutations thereof), corresponding to a large-scale temper- 
ature mode which is sourced by the ISW effect and two small-scale modes which are 
affected by lensing. This closely mimics the local form of non-Gaussianity which is 
predicted in a range of possible early-Universe scenarios (for a review, see [[35]]). For 
upcoming experiments such as Planck, this lensing-induced bispectrum will provide 
a large source of contamination for the measurement of primordial non-Gaussianity, 
and will need to be accurately accounted for to extract the primordial bispectrum pl| , 

Note that because the lensing effects are linear in the observed CMB, lensing 
cannot generate a non-zero bispectrum in the absence of the Ty/ correlation, and the 
lowest order at which the lensing induces non-Gaussianity is the four-point correla- 
tion function or trispectrum. This has been studied thoroughly by several authors [B7L 



5 Lens reconstruction 

In the previous section we discussed how lensing affects the CMB power spectra and 
temperature bispectrum. We turn now to the issue of extracting more information 
about the lensing itself, in particular how to construct an optimal estimator for the 
lensing deflection field from the observed CMB. This reconstruction can then be 
connected to cosmological parameters with a likelihood analysis. 

If we assume that the lensing potential is fixed then the statistics of the lensed 
CMB remain Gaussian, however the inhomogeneities which the lenses introduce 
make it statistically anisotropic. In Fourier space, this manifests itself as correlations 
between modes with 1 ^ I', as we have already discussed. In real space, the manifesta- 
tion is somewhat more intuitive: for an isotropic model, the correlation between two 
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points depends only on the distance between them, and the effect of the lens remap- 
ping is to make the observed, lensed distance between two points slightly larger or 
smaller than the primary, unlensed distance on which the correlation is defined. This 
feature may be exploited to produce estimators for the lensing potential. Consider, 
for example, the lensing expansion in Eq. (|l5|), which shows that at first order in \j/ 
the effect of lensing is to introduce a dependence into the observed CMB on the gra- 
dient of the unlensed CMB. This suggests that a quadratic estimator for the deflection 
angle of the form 

a(x) oc f (x)Vr(x) (34) 

can be used to extract a measure of the lensing signal, if properly normalized. Taking 
the expectation value over realizations of the CMB, this estimator would recover the 
deflection field in the mean (at least to first order in the lensing expansion). The re- 
construction is clearly quite noisy, however, as for any particular realization there will 
be some spurious TVT correlation in the unlensed CMB. Additionally, it is not clear 
how we would obtain the unlensed CMB T(x) with which to perform the correlation. 
It is critical that the gradient of the unlensed CMB be used in Eq. (^) - if instead 
the gradient of the lensed CMB is used, the expectation value over CMB realizations 
vanishes to first order in 

Here we will derive a more rigorous lensing estimator, by maximizing the like- 
lihood of <z(x) given the observed data, loosely following the approach of [|39|]. We 
will take the measurement to consist of noisy maps of T, Q and U measurements 
(or some subset thereof). We combine the observables at each point x into a vector 
@(x) = <9(x) +n(x), where n is the instrumental noise realization, which is also as- 
sumed to be Gaussian. The log-likelihood J?f of the observed CMB is then given by 

-JSf(d|o) = 1 -& T (c M y l + ilndet(C®®), (35) 

where C @0 (x,y) = C @0 (x,y) + N(x,y) is the covariance of the observation (in- 
cluding the effects of lensing) and we have left the integrals over positions implicit. 
We will construct an estimator for the deflection field by maximizing the likelihood 
with respect to the deflection field a(x). Taking the functional derivative of the like- 
lihood with respect to the vector components a, we find 



Sa,(x)-2° {C ' 8a, (x) (C ] ° 2 Tr 



S(^00 



(36) 



The trace term results in a "mean field" over realizations of the unlensed CMB and 
noise. To see this, consider the identity Ti(A) = (x T AC~ l x), where A is any matrix 
and x is a vector of Gaussian random variables with covariance C. Making this sub- 
stitution with C = C 0@ and maximizing the likelihood by setting 5_Sf/5a,-(x) = 
gives the simple equation 



(37) 
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where the components J#f(x) are given by 



&0\-i 







t 8c @e > 
5a,- (x) 



0&\-i 







(38) 



The values of 8C 0@ /8(Xj(x) are calculated analytically from the 0{a) covariance in 
e.g. [ pO[pl[ ] (albeit with a slightly different notational convention than the one used 
here). Equation ([37]) may be solved iteratively using the Newton-Raphson method by 
updating a(x) according to 



5 a, ( j 



dyMj(x,y)[jej-(jrj)}(y), 



(39) 



where sub-scripted quantities are evaluated at the current estimate of a(x), and the 
normalization ^; (x, y) is given by 



/(*,y) = 



-i -i 



5a,-(x)5a / (y) 



5a ; (x) 



(^■-(^■))(y) 



-i 



(40) 



which in practice can be approximated by the Fisher matrix (see below). It is appar- 
ent that the quadratic term Jif(x) produces a measure of the deflection angle a(x). 
Note that in the absence of lensing, with homogeneous sky coverage we must have 
(Jtf(x)) = 0, as the expectation value would otherwise pick out a preferred direction 
on the sky iQ. 

In fact, J^f(x) is quite closely related to the heuristic estimator of Eq. (34). To 
show this, we begin by finding the derivative of the lensed CMB with respect to the 
deflection angle. This is straightforward since 0(x) = 0(x + a(x)), so 



5a; (y) 



6>(x) = 5(x-y)(V,0)(x + a(x)), 



(41) 



i.e. if the deflection angle is changed at a point, the lensed field only changes at the 
same point and by an amount that depends on the gradient of the unlensed fields 
V0 at the source position. To proceed further, it is useful to define the linear lens- 
ing operator A(x,y) such that the lensed CMB is given by <9(x) = 6>(x + a(x)) = 
f c/ 2 yA(x,y)@(y), or explicitly 



A(x,y) = 5(x+a(x)-y). 



(42) 



With this notation the derivative of the lensed covariance 

£,00 = AC ®@ A t (where 

A r (x,y) = A(y,x)) is given by 



^— -C 00 (y,zH5(x-y) J d 2 y>d 2 i! A{y,y')V } ,,C 00 (y',z')A(z,z') + (y <- z ) 



(43) 



Equation (B8J) then simplifies to a real-space operation on filtered fields 0390: 



^(x)=F[(x)(V i F 2 )(x + a(x)), 



(44) 
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where the filtered fields are defined by 

Fi(x) = [(C®®)- 1 6>](x), F 2 (x) = [C®®A T (C®®)- 1 &}(x). (45) 

The filtered field F 2 = A~ l C® @ (C®®y l & is the Wiener reconstruction of the un- 
lensed CMB, i.e. the optimal estimate of the unlensed CMB given the observed noisy 
data, the estimate of a and the unlensed CMB covariance matrix. The gradient term 
(VF2)(x + Oc(x)) is then just the best estimate of the unlensed gradient evaluated at 
the unlensed position, and this is contracted with the inverse-covariance-weighted 
lensed CMB at x to form M\(x). The weighting in F] is typical of any optimal esti- 
mator of CMB observables and serves to down-weight the estimate in areas of high 
instrumental noise or cosmic variance. To understand why the VF2 term is evaluated 
at the displaced position, let the true deflection field differ from the current estimate 
a(x) by 8 a (x). To linear order in 8 a (x), the lensed CMB at x is 

0(x) = 0(x + a(x)) + 8a i (x)(y i 0){x + a(x)) . (46) 

It is therefore necessary to multiply by the gradient of the unlensed CMB at x + a(x) 
to obtain an estimator which is proportional to the desired change in deflection field 
when averaged over realizations of the unlensed CMB. 

We now consider the error on the reconstruction. The Fisher matrix, 



S 2 3? \ 1 



(47) 

provides a lower bound on the variance of an unbiased estimate of a by the Cramer- 
Rao inequality. The expectation value in Eq. ( p7| ) is taken over realizations of the 
CMB and noise for fixed lenses. This may be evaluated analytically at a = if we 
have isotropic sky coverage, in which case we find 

4 0) (x, y ) = 

| £ /V £ /VTr[(C^)- 1 (x',x)V.,C 00 (x,y') (C^)- 1 (y',y)V vJ C 00 (y,x') 



+ (C 00 )- 1 (y,x)V.,C 00 (x, y ') (C®®)- l (y',x')V yJ C®®(x',y) 



■ (48) 

a=0 



This is diagonal in Fourier space in its spatial arguments; for example, for the tem- 
perature alone we have 

(0), v _ f d 2 L 



where N n (L) = - / -, — -~ — ! — TT ! .(49) 

■J K ' 2 J {mf {Cj T +N,) (Cp L| +AT|,_ L |) 

Here, N[ is the noise power spectrum of the temperature map. To make contact 
with previous work jl4] , ^o| ] that reconstructed the lensing potential y rather than 
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the deflection field, we note that the Fisher matrix J^°' is related to that for a by 
i4 0) (x,y) = VyV^-^ 0) (x,y). Equation @ then gives 



#v [x,y) J {my 



1 , rf 2, [L-(L-l)C^ L| +L-lC ; ^i 



2 7 (2^)2 +N l ){C '+N^ M 



(50) 



which agrees with the power spectrum of the reconstruction noise derived for the 
optimal quadratic estimator of y in in the limit of no lensing. Note that in recon- 
structing a rather than y/ we allow for the possibility of a curl mode in a. This can 
act as a monitor for systematic effects, as it is expected to be zero cosmologically. 
However constraints using the iterated estimator can be improved by using a prior 
that it is zero Q. 

We will refer to the Fisher noise level at a = as the "Fisher limit"; it is plotted 
for two instrumental noise levels in Fig. |[ In practice, the reconstruction noise can 
be larger than the Fisher limit due to the confusion created by the lenses themselves, 
which act as a source of noise in the observed CMB and make the lensing likelihood 
slightly non-Gaussian Jl4[]. For temperature and £-mode polarization reconstruction 
these effects are generally negligible, as the lensing represents a small perturbation 
to the primary signal. For a reconstruction including fi-mode polarization measure- 
ments, however, the difference is very important at low noise levels as the lensing 
contributions constitute a significant fraction of the observed fi-mode power. The true 
reconstruction noise can be determined from simulations Jl4]]. In it is shown that 
a good approximation to the true error may be obtained by treating the contribution 
of poorly constrained modes of y/ to be Gaussian. This is done by replacing the un- 



lensed power spectra in the denominator of Eq. (50) with the spectra of the partially 
de-lensed CMB. The amount of de-lensing which is possible is in turn determined 
by the lens reconstruction noise (and so the reconstruction noise calculation must be 
iterated until convergence). We follow this approach to calculate the noise curves in 
Fig. j 

On small scales CMB lensing can also be used to study individual clusters; here 
the unlensed CMB is nearly a pure gradient, and the cluster signal is a small "wiggle" 
in this gradient due to the deflection by the cluster. However the signal-to-noise on 
single clusters is generally low, and for temperature lensing there are several sources 
of confusion (see, for example, Refs. | p3j|44| ]). 



6 Experimental status and prospects 

The experimental evidence for CMB lensing has developed greatly in the past sev- 
eral years. At the power spectrum level, lensing causes ff(l0%) effects at / > 1000 
which are now becoming testable with small-scale experiments. The ACBAR team 
have fit for the amplitude of lensing effects in their parameter analysis of ACBAR and 
WMAP data and find a peak within one-sigma of the expected value (although the 
lensing amplitude is still consistent with zero at the percent level) An optimal 
analysis of lensing effects should produce a more significant detection, but requires 
stringent control of instrumental systematic effects and foregrounds. Thus far, such 
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Fig. 5 Power spectrum of the errors in the gradient part of the reconstructed deflection field for two 
full-sky experiments. Left: an approximation to the Planck satellite, with a(T) = 27 /iK-arcmin. Right: a 
version of the proposed CMBPol mission, with a{T) = 1 pK-arcmin. For both experiments, we assume a 
Gaussian beam of FWHM 7 arcmin and Q and U noise levels equal to \/2a(T). The curves are the theory 
deflection power (solid black), reconstruction errors from temperature alone (red dot-dashed), polarization 
alone (green dashed), temperature and polarization together (black dashed), and the Fisher limit (black 
dotted). For Planck, the Fisher limit is saturated. 



an analysis has only been performed for the temperature data of the WMAP satellite. 
The WMAP data only extends to Z max ~ 500, and so does not contain enough small- 
scale data to reconstruct the lensing potential with good signal-to-noise. By cross- 
correlating a lens reconstruction from the WMAP data with a lower-noise tracer of 
the large-scale structure (LSS), however, it is possible to obtain an overall detection. 
This was first attempted by [[46||, who found no detection of lensing effects in cross- 
correlation with the luminous red galaxies (LRGs) of the Sloan Digital Sky Survey 
(SDSS). Their result was improved upon by [[47|], who used radio sources from the 
NPvAO VLA Sky Survey (NVSS), which has a better redshift overlap with the struc- 
tures which lens the CMB. They found a 3.4c detection of lensing. Finally, [ |48| ] have 
improved their original analysis by incorporating correlation with NVSS sources as 
well as quasars from the SDSS survey. They use a more conservative cut of the NVSS 
sources and a slightly sub-optimal estimator compared to j47[], obtaining evidence for 
the expected lensing-LSS correlation of 2.5c. All three of these cross-correlation re- 
sults have incorporated very thorough tests for both instrumental and astrophysical 
systematics, most of which have negligible effects on the current overall error budget. 

At present, therefore, there is experimental evidence at modest significance for 
the effect of gravitational lensing in the CMB. High-resolution experiments which 
are currently taking data, such as Planck, ACT, and SPT, are expected to improve 
this situation dramatically, with reconstruction signal-to-noise approaching unity for 
many modes of the \j/ field. This will enable CMB lensing to become a useful tool 
for precision cosmology. With lens reconstruction from the Planck data, for example, 
constraints on the energy density of massive neutrinos are expected to improve by a 
factor of two over an analysis which neglects a lensing reconstruction ||J|(j.[t9|]. CMB 
lensing measurements will also enable definitive testing between some theoretical 
models of the anomalous CMB "cold spot" J50|]. Upcoming large-scale polarization 
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experiments, such as BICEP2, QUIET, and PolarBeaR, are also expected to detect 
lensing in polarization through the B-modes which it induces. 

As CMB lensing moves from a purely theoretical study to an observational one, 
there are a number of challenges to be overcome. The effects of foregrounds, such as 
the Sunyaev-Zel'dovich effect and residual extra-Galactic point-sources, will become 
increasingly important as high signal-to-noise detections of lensing are attempted. 
The construction of a realistic lensing likelihood for the purpose of parameter con- 
straints has also yet to be demonstrated. In the more distant future, CMB lensing 
should provide the dominant source of confusion for primordial Z?-mode measure- 
ments targeting the signal from recombination around / ~ 100 (for sufficiently low 
tensor amplitude). Lens reconstruction with polarization and delensing will become 
important tools to mitigate this confusion, and preliminary work on these techniques 
and their sensitivity to foregrounds and instrumental systematics is currently under 
active development . 
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